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Generating entanglement by simply cooling a system into a stationary state which is highly entan- 
gled has many advantages. Schemes based on this idea are robust against parameter fluctuations, 
tolerate relatively large spontaneous decay rates, and achieve high fidelities independent of their 
initial state. A possible implementation of this idea in atom-cavity systems has recently been pro- 
posed by Kastoryano et al. [Phys. Rev. Lett. 106, 090502 (2011)]. Here we propose an improved 
entanglement cooling scheme for two atoms inside an optical cavity which achieves higher fidelities 
for comparable single- atom cooperativity parameters C . For example, we predict fidelities above 
90% even for C as low as 20 without requiring individual laser addressing and without having to 
detect photons. 



PACS numbers: 03.67.-a, 42.50.Pq 



I. INTRODUCTION 

Current atom-cavity experiments with coupling con- 
stants cavity decay rates and atomic decay rates 
F operate in a parameter regime where the single-atom 
cooperativity parameter C, 



C 



(1) 



is at most one or two orders of magnitude larger than one 
[IHQ. However, the practical realisation of atom-cavity 
quantum computing schemes usually requires Cs above 
200 to achieve single-operation fidelities above 90 % [7- 
[T4] . The only alternative are probabilistic quantum com- 
puting schemes. These promise fidelities above 90 % even 
when C = 10 but rely either on the detection of single 
photons [15l |16] or on the observation of macroscopic flu- 
orescence signals [TT. Because of being conditional, they 
require relatively high photon detection efficiencies and 
cavity mirrors with low absorption coefficients. Using 
currently available experimental setups to entangle atoms 
in optical cavities with a very high fidelity therefore re- 
quires a different approach than previously proposed in 
the literature. 

Recently it has been pointed out by several authors 
p!8W26| that it is possible to generate entanglement in a 
controlled way by simply cooling qubits into well-defined, 
highly entangled states. The main idea behind this ap- 
proach is to design laser fields such that the target state 
becomes the stationary state of the system. State prepa- 
ration schemes based on this idea are expected to toler- 
ate much higher spontaneous decay rates than proposals 
which do not use dissipation in this way. Moreover, when 
cooling a system into an entangled state, the fidelity of 
the state preparation no longer depends on the initial 
state of the system which makes the entanglement gen- 
eration more robust against errors. Although being very 
promising, this approach has only recently been studied 
as a tool to entangle two atoms trapped inside the same 
optical cavity. The only examples are Kastoryano et al. 




FIG. 1: (a) Experimental setup to cool two atoms inside an 
optical cavity into a maximally entangled state. Here, F and 
n denote the spontaneous atom and cavity decay rates while 
Qo and Qi are the relevant laser Rabi frequencies, (b) Level 
scheme of a single-atom. The 1-2 transition couples reso- 
nantly with coupling constant g to the cavity field. The spon- 
taneous decay rates for the 0-2 and the 1-2 transitions are 
Fo and Fi with T = Fo + Fi. 



[27] and Wang and Schirmer [28^. 

In this paper we follow similar ideas as Refs. [271 [M] 
and design an entangling scheme to cool two atoms in- 
side an optical cavity into a maximally entangled state. 
As proposed in Ref. [23 , and in close analogy to the 
laser sideband cooling technique of trapped ions [29 , we 
employ level shifts and apply laser fields such that only 
the target state experiences off-resonant driving. Every 
ground state of the system other than the target state 
couples resonantly and sufficiently strong to rapidly de- 
caying excited states. Doing so, the target state becomes 
the stationary state of the quantum system. It is reached 
independently of the initial state of the system after a 
certain transition time. As in laser sideband cooling, the 
fidelity of the final state reaches one when the detuning 
of the target state becomes much larger than the relevant 
laser Rabi frequencies and decay rates. 

The concrete experimental setup which we consider in 
this paper is shown in Fig. [T] It consists of two atoms si- 
multaneously trapped inside an optical cavity. The main 
decay channels in this system are the spontaneous emis- 



2 



sion of photons from the excited atomic state |2) with 
decay rate F and the leakage of photons through the cav- 
ity mirrors with decay rate k,. Both atoms are driven 
by external laser fields which couple respectively to the 
0-1 and the 1-2 transition. In the following, we design 
the detunings and Rabi frequencies of these laser fields 
such that the stationary state of the atom-cavity system 
is to a very good approximation given by the maximally 
entangled atomic ground state 



1+) = (|01)+|10))/^/2, 



(2) 



while there is no photon inside the cavity. As we shall 
see below, individual laser addressing of the atoms is not 
required. The entangling scheme proposed in this paper 
uses an energy shift of the target state which is due to 
a non-zero atom-cavity coupling constant g as well as 
spontaneous emission from excited states. This makes it 
possible to prepare the state in Eq. Q with a fidelity 
above 90% even when C is as low as 20 and without 
having to detect photons. 

One advantage of the state preparation scheme pre- 
sented in this paper is that it predicts higher fideli- 
ties than the recently proposed entangling schemes in 
Refs. [271 |28 although they employ similar level shifts 
to cool two the atoms into a maximally entangled state. 
Ref. [2T uses a similar atomic level scheme as the one 
shown in Fig. [T] but with the addition of a driven mi- 
crowave transition between the triplet states. Ref. [28 
relies on the presence of an external magnetic field gra- 
dient to produce the required level splittings. 

There are five sections in this paper. In the next sec- 
tion, we introduce a four-level toy-model and discuss how 
to cool it into one of its ground states. The reason for 
the introduction of this toy-model is that the entangling 
scheme proposed in this paper cannot be modelled eas- 
ily analytically. There is no interaction picture in which 
the system Hamiltonian becomes time-independent. Al- 
though being much simpler, the toy-model in Section 
\n\ captures all the basic features of the proposed state 
preparation scheme, provides much insight into its cool- 
ing mechanism, but is nevertheless analytically tractable. 
In Section [nil we present all the details of our entangling 
scheme, draw analogies to the toy-model, and support 
our claims about the parameter dependence of its fidelity 
with the help of numerical simulations. We finally sum- 



marize our findings in Section IV 



II. STATE PREPARATION IN A TOY-MODEL 

In this section we consider a simple four-level system 
and pose the task to prepare it in one of its two ground 
states. The role of the experimental parameters in this 
simple model, i.e. its laser Rabi frequencies, detunings, 
and spontaneous decay rates, can later be mapped onto 
the atom-cavity coupling constant the cavity and the 
atom decay rates and F, and laser parameters Qi and 




FIG. 2: The toy-model level scheme. The 0-2 and the 1-3 
transitions are driven by a laser field with Rabi frequency Q 
and a detuning A with respect to the 0-2 transition. The 
excited atomic states both decay spontaneously into |0) and 
|1) with a decay rate F/2. 



understanding of the toy-model will allow us to correctly 
predict the general dependence of the fidelity and the 
cooling rate of the proposed entangling scheme on these 
experimental parameters. 



A. Theoretical model 

The toy- model contains only two ground states |0) and 
|1) and two excited states |2) and |3), as shown in Fig. [5] 
For simplicity, we assume that the decay rates for all 
four decay channels are the same and denote the overall 
spontaneous decay rate of level 2 and 3 by F. Moreover, 
we assume that the system is driven by a single laser 
field of frequency cjl and Rabi frequency Q. This laser 
is in resonance with the 1-3 transition but detuned from 
the 0-2 transition by a detuning A. The spontaneous 
emission of photons is in the following taken into account 
using the master equation 



(3) 



Si in the entangling scheme proposed in Section III Our 



where g is the density matrix of the system. The con- 
ditional Hamiltonian i^cond describes the time evolution 
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under the condition of no photon emission, while R{g), 

^(^) = E E l^\m9\m, (4) 

i=0,l j=2,3 

relates to the reset state in case of a photon emission. 
Within the rotating wave approximation and in the in- 
teraction picture with respect to the free Hamiltonian 



Ho = ^tuvi\i){i\-hA\2){2\ , 



(5) 



i=0 



where huJi is the energy of level i^cond equals 
ffcond = ^/l!^(|0X2| + |l)(3|+H.c.) + fiA|2)(2| 
-V(|2)(2| + |3)(3|), 



(6) 



since cjl = 0J2 — A. 

The master equation in Eq. (|3| can now be used to 
calculate the fidelity of the proposed state preparation 
scheme, i.e. the stationary state population of its tar- 
get state |0), and its cooling rate. This is most easily 
done using rate equations which are a complete set of 
differential equations for the time evolution of expecta- 
tion values. The time derivative of an expectation value 
of a time-independent operator A equals 

(A) = Tv{Ag). (7) 

The above master equation hence implies that 

(A) = -iif^{[A|0){2| + |l)(3|+H.c.]) 

-iA{[A,|2)(2|])- ^ lr{A\j){j\ + \j){j\A) 

i=2,3 

+ E E lr{\j){i\A\i){j\). (8) 

i=0,l j=2,3 

In the following we consider the Hermitian operators 
K)01 + b')(^l^ and - and denote their 

(real) expectation values by 

Pi = {i\Q\i). 

kij = 2Im(z|^|j), 

kj = 2Re {i\g\j). (9) 

Substituting these operators into Eq. (Is]) yields 



Po 

Pi 



-^nko2 + ^r{P2 + P3) , 

lnhs + ^T{P2 + P3) , 



1 



P2 = -nko2-rP2, 
P3 = ^nkis-rPs, 
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^r/co2 . 



13 ■ 



(10) 



These seven equations form a complete set of rate equa- 
tions and are sufficient to analyse the time evolution of 
the population Pq in the target state |0) which equals the 
fidelity F. 



B. The basic idea 

Suppose we aim to transfer the toy- model in Fig.[2]into 
one of its ground states, for example the |0) state, with 
a very high fidelity and without having to control the 
initial state of the system. This is possible when laser 
driving is applied such that the |0) state becomes the 
stationary state of the toy-model as it applies when the 
laser detuning for the 0-2 transition is much larger than 
the other system parameters, i.e. when 



A > 1], r. 



(11) 



This condition guarantees that it is much more likely 
for the system to spontaneously decay into the |0) state 
when being in one of the other three states, than being 
driven out of it [23]. Once the system has reached its 
stationary state it therefore remains there with a very 
high probability. 

This is confirmed by Fig. |3] which shows the time de- 
pendence of 1 — F, i.e. of the total population in states 
other than the target state |0), for a wide range of ex- 
perimental parameters. The solid lines in this Fig. |3] are 
the result of a numerical integration of the rate equa- 
tions in Eq. (10) which assumes the worst case scenario 



with the toy- model initially in The plots show expo- 
nential cooling towards the target state until the system 
reaches a stationary state. The fidelity of the state prepa- 
ration equals the population of the |0) state and is indeed 
very close to unity, as long as condition ( 11 ) applies. The 



cooling rate and the fidelity of the state preparation both 
depend on the relative size of Q and F with respect to 
the detuning A. 



C. Stationary state fidelity 

To identify the best way of preparing the target state, 
we now derive approximate analytical expressions for the 
stationary state fidelity, F, and the cooling rate, 7c, which 
is a measure for the time it takes the system to reach its 
stationary state. Since F is the stationary state popula- 
tion Po in the |0) state, it can be calculated simply by 
setting the time derivatives of the expectation values in 
Eq. (10) equal to zero. Doing so, we find that 



3n^ 



4A2 + 4^]2 + 



(12) 
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FIG. 3: Logarithmic plot of the time dependence of the dis- 
tance 1 — F from the target state |0) for different values oiQ/ A 
and r/A. The system is initially in In the upper plot we 
have r = 0.2 A. In the lower plot we have Q = 0.05 A. The 
solid lines are the numerical solutions of the rate equations in 
Eq. ( |10|). The dashed lines illustrate the analytical solution 
in Eq!"^T]). 



For relatively large detunings A, as in Eq. ( 11 ), this equa- 
tion simplifies to 



1 



4A2 



(13) 



This result confirms that the fidelity is close to one in the 
parameter regime given by Eq. ( 11 ). Fig. [4] illustrates the 
effects of finite ft and F. It also shows that an increase in 
the Rabi frequency Q reduces the stationary state fidelity 
more rapidly than an increase in the decay rate F. 



D. Heating and cooling rates 

To see how quickly the toy-model reaches its station- 
ary state, we now introduce the notion of a cooling and 
a heating rate which we denote 7c and 7h, respectively. 
For simplicity, and since we are anyway only interested in 
the general scaling of these rates with the experimental 
parameters, we assume that these rates do not depend 
on the current state g of the system. Invoking the con- 
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FIG. 4: Contour plot which shows the stationary state fidelity 
F in Eq. Q as a function of 0/A and F/A. 



servation of probability flux, we then find that 
Po = lc (l-Po)-7hPo. 



(14) 



The principle here is that the rate at which the fidelity, 
i.e. the current population in the |0) state, changes in 
time is equal to the rate at which population is cooled 
into the target state minus the rate at which population is 
heated out of the target state. When the system reaches 
its stationary state, the fidelity remains constant. The 
above equation hence implies 



7hF = 7c (1-F) . 



(15) 



Since we already know F (cf. Eq. (12)), this relation can 
be used to obtain the cooling rate after obtaining an es- 
timate for the heating rate. As we shall see below, it is 
easier to derive an approximate expression for 7h, than 
calculating 7c directly. 



Considering the parameter regime in Eq. ( 11 ), the rate 
equations in Eq. (10) can be simplified via an adiabatic 
elimination. Only the coherences ko2 and /o2 evolve on 
the fast time scale given by A. Setting their time deriva- 
tives equal to zero, we find that 



2Tn 



^02 



r2 + 4A2 



(16) 



Assuming that the toy- model is in |0), i.e. that Pq = 1 
and Pi = P2 = P3 =0, and substituting the above 
expression for /co2 into the rate equation for Pq yields 



Po = 



r2 + 



(17) 
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FIG. 5: Contour plot which shows the cooling rate 7c in 
Eq. ([l9l as a function of Q/A and T/A. 



Comparing this equation with Eq. (14) for Pq = 1, we 



find that the heating rate is to a very good approximation 
given by 



7h 



r2 + 4A2 



(18) 



This is the rate at which the target state |0) loses its 
population. Substituting this result and Eq. (12) into 
Eq. (|T5|, we get 



(4A2 + r2)(3fi2 + r2) 



(19) 



Fig.|5]shows this coohng rate 7c for a wide range of exper- 
imental parameters. For relatively small Rabi frequencies 
the cooling process becomes faster with increasing Q. 
However, it is not worth increasing Q beyond a certain 
size which saturates 7c. In the parameter regime given 
by Eq. (11), the cooling rate 7c simplifies to 



7c 



(20) 



which no longer depends on A but only holds for suffi- 
ciently large detunings. 

In order to get a feeling for the accuracy of the cooling 
rate 7c in Eq. (19), we now solve Eq. (14) analytically 



and compare the result with exact numerical solutions of 



the rate equations in Eq. (10). Doing so and assuming 
Po(0) = we find that 



7c 



7c + 7h 



1 



-(7c+7h)t 



(21) 
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FIG. 6: (a) Logarithmic plot of the time dependence of the 
distance 1 — F from the target state |0) for the case where the 
system is initially in |1) and where F = 0.2 A and Q — 0.05 A. 
(b) Same as (a) but for a time dependent laser pulse with a 
Rabi frequency Q{t) as in Eq. (22) with Qq = 0.05 A. (c) 
Theoretical minimum for 1 — F obtained from Eq. (23) for 
F = 0.2 A. 



Fig. |3] compares this analytical result with numerical 
solutions of Po(^) for different experimental parameters 
Q/A and F/A. It shows that Eq. (19) reflects the gen- 
eral parameter dependence of the cooling rate on Q/A 
and F/A correctly. The above approximate solution is in 
general slightly higher than the actual cooling rate. The 
reason for this is that the heating rate in Eq. (18) has 



been calculated for the case, where the system is initially 
in |0), i.e. when it is the highest. 



E. Choosing experimental parameters 

Fig. [4] shows that maximising the stationary state fi- 
delity F requires a Rabi frequency Q as small as possible. 
However, from Fig. [5] we see that we only obtain high 
cooling rates when Q is relatively large. To minimise the 
state preparation time while maintaining a high fidelity, 
we therefore suggest using a laser pulse with a time- 
dependent Rabi frequency to prepare the target state. 
This laser pulse should be large initially and should reach 
zero by the end of the cooling process. For example one 
could choose 



3^0 



(l + 7c(0)t)^ 



(22) 



with 7c (0) being the cooling rate in Eq. ( |2Q[ ) for the ini- 
tial Rabi frequency Q{0) = Qq. Alternatively, one could 
choose an exponentially decreasing Rabi frequency. How- 
ever, in this case, Q would drop off too rapidly, thereby 
resulting in a fidelity that is far from optimal. Here we 
do not discuss how to optimise the spontaneous decay 
rate, since F is in general fixed. 

Fig. [6| co nfirms that choosing the Rabi frequency Q as 
in Eq. (22) indeed yields a significant speed up compared 
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to time-independent Rabi frequencies. We also observe a 
stationary state fidelity which is close to the theoretical 
maximum obtained when setting = 0. From Eq. (13) 



we see that this maximum is to a very good approxima- 
tion given by 



F{n = 0) = 1 



4A2 + 2r2 • 



(23) 



It is indicated by a dashed line in Fig. [6| In the next Sec- 
tion we use similar time-dependent laser pulses to prepare 
two atoms inside an optical cavity relatively fast and with 
a high fidelity in a maximally entangled state. 



III. ENTANGLING SCHEME 

The toy-model described in the previous section is 
based on a simple principle: driving populations out of 
all undesired states resonantly while driving the target 
state off-resonantly. This approach can indeed be used 
to prepare target states with a high fidelity [23 . In this 
section, we use this idea to prepare two atoms inside 
an optical cavity (cf. Fig. [T]) in the maximally entangled 
state 1+) in Eq. The first half of this section presents 
a theoretical description of the atom-cavity system. Af- 
ter identifying its dressed states, we select appropriate 
laser Rabi frequencies and detunings. As we shall see be- 
low, the state preparation requires three different driving 
lasers but there is no need to address atoms individually. 
A comparison with the toy-model introduced in the pre- 
vious section allows us to predict the dependence of fi- 
delity and cooling rate of the proposed state preparation 
scheme on the experimental parameters with a very high 
accuracy. 



A. System Hamiltonian without laser driving 

The experimental setup which we consider in this pa- 
per consists of two atoms placed inside an optical cavity 
as shown in Fig. [l] Each atom contains a A-type level 
configuration with huj and \j) denoting the correspond- 
ing energies and energy eigenstates {j = 0, 1, 2). Suppose 
the 1-2 transition of each atom couples resonantly with 
coupling strength g to the quantised cavity field mode 
with frequency cj, 
system equals 



Then the Hamiltonian i^gys of this 



2 2 2 

H,y, = ^%|l),,(2|ct+H.c. + ^^fz£^, 

i=l i=l j=0 

-^huJcC^c (24) 

in the absence of external laser driving. Here c and 
are the cavity photon annihilation and creation operators 
for a single photon inside the optical cavity. As we shall 
see below, it is important that the atomic states |0) and 



Energy eigenstate 


Energy 


|00,0) 





|+,0)=(|01,0)+|10,0))/V2 




|-,0)=(|01,0)-|10,0))/V2 




|11,0) 


2nwi 



TABLE I: Energy eigenstates a nd energy eigenvalues of the 
system Hamiltonian i^gys in Eq. p4| for the Ho subspace with 
no atom in |2) and no photons in the cavity. 



|1) differ in energy by an amount which is significantly 
larger than hg. 

In the following, we identify the relevant energy eigen- 
states of this Hamiltonian, since this will allow us to iden- 
tify appropriate laser drivings and detunings for the pro- 
posed state preparation scheme. To do so, we denote 
states with atom 1 in atom 2 in \j2) and n photons 
in the cavity by \jij2^n)- Moreover we notice that the 
Hamiltonian i^gyg preserves the total amount of popula- 
tion in the excited atomic state |2) and the cavity field 
mode. It therefore acts on fixed excitation subspaces Hn 
of the complete Hilbert space whose energy eigenstates 
can be calculated separately. The eigenstates and eigen- 
values of the subspace Hq of states with no population of 
the excited atomic state |2) and no photons in the cav- 
ity are summarised in Table |Tj Table |TT] shows the eight 
energy eigenstates and the corresponding energy eigen- 
values of the subspace Hi of states with either one atom 
in 1 2) or one photon in the cavity and adopts the notation 



iMi) 
IMo,±) 
|Ao,±) 

|Ai,±) 



(|21,0)-|12,0))/V2, 

(|02,0)-|20,0)±|01,1)t|10,1))/2, 

(|02,0)+|20,0)±|01,l)±|10,l))/2, 

(|12,0)+ |21,0)± V2|ll,l)) /2. (25) 



Fortunately, there is no need to identify the energy eigen- 
states of the atom-cavity system of the subspace of states 
with more than one excitation in the atomic state |2) and 
the cavity field mode. The reason for this is that these 
states do not couple directly to the states in Table |l] in 
case of laser driving. Therefore they do not have to be 
taken into account when choosing laser parameters such 
that only the |+,0) state experiences off-resonant laser 
driving. 



B. Laser driving 

As we shall see below, the state preparation of the 
maximally entangled atomic state |+) in Eq. Q requires 
the simultaneous excitation of the two atoms with three 
different laser fields. In the following we assume that 
the 0-2 transition of each atom is driven by two differ- 
ent lasers with Rabi frequencies Qq^^ and frequencies uJq^^ 
respectively {k = 1,2). The 1-2 transition of each atom 
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Energy eigenstate 


Energy 


|00,1) 


huc — h{uj2 — CJi) 


iMi) 


h{(jUi + CJ2) 


|Mo,±) 


h{uj2 ± 9) 


|Ao,±) 


h{(jj2 ± g) 


|Ai,±) 





TABLE II: Energy eigenstates a nd e nergy eigenvalues of the 
system Hamiltonian ifsys in Eq. ( p4| for the 1-Li subspace of 
states with either one atom in |2) or one photo n in the cavity. 
The table uses the notation introduced in Eq. (25). 



should moreover be driven by a laser field with Rabi fre- 
quency r^i and frequency cjli. The laser Hamiltonian 
in the Schrodinger picture and the usual rotating wave 
approximation is then given by 



2 2 



i=i fc=i 

+ ^-;ii^ie"^"*|l),i(2|+H.c. (26) 



The realisation of this Hamiltonian does not require in- 
dividual laser addressing, since both atoms experience 
exactly the same laser driving. 

In order to see how to best choose the laser frequen- 
cies co^Lo ^Li, we now consider the effect of this laser 
Hamiltonian on the Ho subspace. This effect can be de- 
scribed by the restricted laser Hamiltonian Hi^ defined 
as 



(27) 



with the projector P being the projector on Hq and Hi 
given by 



(28) 



Using the eigenvectors of the undriven atom-cavity sys- 
tem Hamiltonian which can be found in Tables [J and [III 
one can show that this Hamiltonian equals 



|00,0)(Ao,x| 



2 ^ 

k=l 

E -^fil^ie'-"*[|+,OXAo,x| 



2V2 

+ \-,0){no,x\ + |ll,OXAi,a; 



Ground 
state 


Excited 
state 


Rabi 
frequency 


Effective 
detuning 


|oo,o> 


|Ao,±) 
|Ao,±) 




(1) _i_ 

(2) 1 


l+,0) 


|Ai,±> 
|Ai,±) 
|Ao,±) 




^Lo - ^2 ± V2g 
^ho - ^2 ± V2g 


l-,o) 


iMi) 

Imi> 

Imo,±> 


Qi 


(1) 

^LO - ^2 

(2) 

^LO - ^2 

c^Li + (^1 - (^2 T g 


|11,0> 


|Ai,±> 


Qi 


c^Li + cji - T V^g 



TABLE III: Most relevant laser-driven transitions of the 
atom-cavity system in the dressed state picture. The table 
shows the respective ground and excited states and indicates 
the corresponding laser parameters. 



possible, since there are more laser fields than atomic 
transitions. 

To make it nevertheless easy to identify the relevant 
laser Rabi frequencies and detunings, we now transform 
the laser Hamiltonian ^l(^) for the subspace Ho0Hi into 
the interaction picture with respect to the system Hamil- 
tonian in Eq. ( 24 ) . Taking into account the eigenvalues of 
this Hamiltonian which can be found in Tables HI and Ull we 
obtain another time-dependent Hamiltonian from which 
we can directly read off the information which is relevant 
for the construction of an entangling scheme via cooling. 
The result of this calculation is summarised in Table Hill 
which shows all laser-driven transitions and states the 
corresponding relevant laser parameters. 



C. Effect of spontaneous emission 

As has been illustrated already in Section [TTj dissipa- 
tion is an essential component of state preparation via 
cooling. In the atom-cavity system analysed in this sec- 
tion, dissipation can occur via the photon emission from 
the excited atomic state |2) with the spontaneous decay 
rate F and via the leakage of a photon through the cavity 
mirrors with the spontaneous decay rate k,. The condi- 
tional Hamiltonian that describes the time evolution of 
the atom-cavity system between photon emissions equals 

2 

Hcond = H,y,^Hi^{t)-^hT^\2)u{2\-^hKc^c. 

(30) 

The first two terms in this equation are the system Hamil- 



tonian Hsys in Eq- (24) and the laser Hamiltonian Hi,(t) 



in Eq. (26). In case ofan emission, the density matrix of 



H.c. (29) the atom-cavity system changes up to normalisation into 



in the Schrodinger picture. Changing into an interaction 
picture in which Hi,(t) becomes time independent is not 



^(^') = E E rjli)«(2|p|2)i,01 + Kcpct, (31) 



j=0,l i=l,2 
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Ground 
state 


Excited 
state 


Rabi 
frequency 


Effective 
detuning 


|00,0) 


|Ao,+> 
|Ao,-) 
|Ao,±) 




-25 



l+,o> 


|Ai,±) 
|Ai,±> 




±(72±1)3 




|Ao,±) 




-{V2±l)g 


|-,o> 


Imi> 

iMi) 




-9 





Imo,±> 




-(V2±l)5 


|11,0) 


|Ai,+> 


Qi 






|Ai,-> 


Qi 






TABLE IV: Transitions between dressed states driven near 
resonance by the application of three lasers with Rabi fre- 
quency ^^o"'^'*, r^o^'* and Qi. 



where Tj denotes the spontaneous decay rate of the 
atomic 2-j transition. The overall decay rate of the ex- 
cited atomic state is given by T = Fq + Fi. Overall, the 
time evolution of the system in the presence of sponta- 
neous emission is described by master equations which 
are of exactly the same form as the master equations in 
Eq. (§. 

D. Appropriate laser parameters 

As already mentioned above, the target state of the 
state preparation which we propose here is the maxi- 
mally entangled atomic state |+) in Eq. ([2|. In order to 
assure that this state becomes the stationary state of the 
atom-cavity system in Fig. [l] we need to choose the laser 
frequencies oo^q ^ ^lo^ such that the |+,0) expe- 

riences only off-resonant driving, while the states |00,0), 
| — , 0), and 1 11, 0) couple resonantly to at least one of the 
three driving lasers. Having a closer look at Table |III[ 
we see that this applies, if we choose 

(1) 

^LO = ^2 - ^ , 
(2) 

^LO - ^2 , 

^Li = 002 — (^1 — V2g . (32) 

Table |IV| shows the effect of this choice of laser frequen- 
cies on the sixteen transitions which need to be taken into 
account when designing the state preparation scheme 
proposed in this paper. 

The system Hamiltonian H^ys treats both atoms in ex- 
actly the same way. Its eigenvectors are therefore ei- 
ther symmetric or antisymmetric with respect to an ex- 
change of the two atoms. The same applies to the ef- 
fective laser Hamiltonian Hh{t). Since both atoms ex- 
perience exactly the same Rabi frequencies, the lasers 
excite either transitions between two symmetric states 



|Ai,+> 



y/2g 




FIG. 7: Level configuration showing the laser driving, Rabi 
frequencies, and detunings experienced by the target state 
|+,0) in the dressed state picture. For simplicity, we show 
only the least detuned couplings. 



or two anti-symmetric states. This allows us to consider 
the symmetric and the antisymmetric state space sepa- 
rately when analysing the effect of the laser driving in the 
dressed state picture of the atom-cavity system. There 
are three symmetric ground states and one antisymmet- 
ric ground state. These are {|00,0), |+,0), |11,0)} and 
{| — ,0)} respectively. 

Fig. [7| illustrates the laser driving experienced by the 
target state |+,0). Since this state is a symmetric state, 
the relevant level configuration involves only the target 
state and the four symmetric states with one excitation 
in 1 2) or the cavity mode. As one can see from Table 
[TVl in the dressed state picture, these three lasers in- 
volve |+, 0) in six different transitions. For simplicity, we 
show only the least detuned couplings for each laser. We 
see that the target state |+,0) experiences indeed only 
off-resonant driving. The smallest and therefore most 
relevant detuning is given by 

Synin = iV2-l)g, (33) 

as long as the frequency uJi is sufficiently larger than the 
atom-cavity coupling constant g. All other states with no 
excitation are resonantly driven by one laser field. This 
is illustrated in Fig. [S] which shows the resonant transi- 
tions in the symmetric and the antisymmetric subspace 
separately. One laser couples 1 00, 0) to |Ao, — ) with Rabi 
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FIG. 8: Level configuration siiowing all the resonantly driven transitions in the dressed state picture, their Rabi frequencies, 
and their detunings in the subspace with zero or one excitation in |2) or the cavity mode, (a) Symmetric subspace with 
ground states |00,0), | + ,0), and |11,0). (b) Antisymmetric subspace with ground state | — ,0). 



frequency q!q^/^/2. Another laser couples | — ,0)to 

with Rabi frequency q!^^ /^/2^ while a third laser drives 
|11,0) into |Ai, — ) with Rabi frequency Vti. In principle, 
we would like these lasers which empty unwanted states 
to be relatively strong. However, it is not possible to in- 
crease them without increasing also the Rabi frequencies 
for the off-resonant driving of the target state shown in 
Fig. [71 



E. 



Fidelities and cooling rates for constant laser 
driving 



Eq. (33) shows that the minimum detuning experi- 
enced by the target state (^min depends only on the atom- 
cavity coupling constant g. Eq. (11) in Section |ll| there- 
fore suggests that the stationary state of the atom-cavity 
system in Fig. [l] is to a very good approximation given 
by the state |+,0) as long as 



(1) o(2) 



cJi > ^ > F, A., % \ Vl, 



' 



(34) 



In other words, for this parameter regime we can expect 
the atoms to be with a very high fidelity in the maximally 
entangled state |+) in Eq. ([2| after a certain transition 
time t. A comparison with the toy-model state prepa- 



ration scheme in Section |TI| even yields approximate so- 
lutions for the fidelity F and the cooling rate 7c of the 
proposed entangling scheme. 

For simplicity, we assume in the following that all 
three cooling lasers have the same Rabi frequency 
and that the two atomic decay rates, Fq and Fi, are 
equal. The only remaining spontaneous decay rates are 
the spontaneous atom decay rate F and the cavity pho- 
ton leakage rate k. For example, the symmetric state 
(|Ao, +)+|Ao, —)) /a/2 with one atom in |2) has the sponta- 
neous decay rate F, whilst the state (|Ao, +)— |Ao, — ))/a/2 
with one photon in the cavity decays with k. We in- 
fer from this that the fidelity of the proposed entangling 
scheme depends on the size of both decay rates. Taking 
this into account, we replace the spontaneous decay rate 
F in Eq. (20) in the following with the average of k and 

r, 



1 



(A. + F). 



(35) 



The analog of the laser detuning A in the toy-model state 
preparation scheme is the detuning (5min in Eq. ( 33 ) which 



is the minimum laser detuning experienced by the target 
state |+,0) during the cooling process. Taking into ac- 
count that 

A (^^in, (36) 
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1.0 r 




1.0 r 



FIG. 9: Stationary state fidelity F of the proposed entangling 
scheme for different spontaneous decay rates n and F and 
Q = 0.03^ obtained from a quantum jump simulation. Here 
K and F are varied such that the cooperativity parameter C 
remains constant at C = 25. 




FIG. 10: Stationary state fidelity F of the proposed entangling 
scheme as a function of the cooperativity parameter C for 
Q = 0.03^ and n — 2F. (a) Numerical solution of the time 
evolution of the system using the quantum jump approach, 
(b) Analytical result in Eq. (37). 



Eqs. (20) and (13) suggest that 7c and F are to a very 



good approximation given by 

2i]2(^ + r) 



7c 



F = 1 



16(V^-l)V 



(37) 



This result confirms Eq. ( |34[ ) which suggests that fideli- 
ties F close to one are only obtained when g is much 
larger than all other system parameters. 

The remainder of this paper confirms these approxi- 
mate solutions with the help of a numerical analysis of 
the proposed entangling scheme. The following analysis 
is based on the quantum jump approach [30 which allows 
us to simulate all the possible trajectories of the atom- 
cavity system in Fig. [T] By averaging over many trajec- 
tories, we obtain an approximate solution of the master 
equation in Eq. ([3| . To calculate the no-photon time evo- 
lution of the system we use the conditional Hamiltonian 
i^cond in Eq. ([30|. In case of a photon emission we reset 
the atom-cavity system such that its state after a photon 
emission is on average given by 7l{g)At with 7l{g) as in 
Eq. (31). For simplicity, we consider only relatively small 



Rabi frequencies. In this way we avoid the population of 
highly excited states. The population of such states is 
not expected to decrease the fidelity of the final state 
since they decay relatively rapidly. However, in this way 
we can restrict the size of the Hilbert space which has to 
be taken into account during simulations to states with 
at most three photons in the cavity. 

Let us first have a closer look at how changing the 
relative size of hi with respect to F affects the fidelity 
F of the state preparation. Fig. |9] shows F for differ- 
ent spontaneous decay rates hz and F. These are chosen 
such that the single-atom cooperativity parameter C in 
Eq. ([T]) remains constant at C = 25. Each data point 



represents the average fidelity calculated from a time se- 
ries like the one shown in Fig. [TT] As a result we find 
that the proposed state preparation scheme works best 
when T — = 0.15 ^. This implies that ideally one should 
have 



2r, 



(38) 



when C = 25. We therefore assume that this applies in 
the remainder of this section. 

The main result of this subsection is an estimation of 
the stationary state fidelity of the maximally entangled 
atomic state |+) which can be achieved with the pro- 
posed entangling scheme. To establish this numerically, 
we use a series of time evolutions such as those shown 



in Fig. 11 and average over the fidelity once the system 
is approximately in its stationary state. Fig. 10 shows 
the stationary state fidelity of the target state as a func- 
tion of the cooperativity parameter C for a constant laser 
Rabi frequency Q. The numerical results are compared 
with the analytical result for F in Eq. (37). Indeed we 



find very good agreement between analytical and numer- 
ical results. It is clear from Fig. [lO] that the achievable 
fidelity F increases rapidly with increasing cooperativity 
parameter C. However, fidelities above 90% are possible, 
even for a cooperativity parameter C as low as 20. 



F. Minimising the state preparation time 

As in Section [TTj we find that a relatively large cooling 
rate 7c requires relatively large Rabi frequencies. At the 
same time, we only obtain a fidelity F close to unity 
for very small Rabi frequencies. In order to maximise 
the fidelity of the state preparation while maintaining 
a substantial cooling rate, we therefore proceed in the 
following as in Section [II E| and assume a time dependent 
Rabi frequency Q. Similarly as in Eq. (22), we assume in 
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(a) 
-(b) 



200 400 600 800 1000 

t in units of 1/g 

FIG. 11: (a) Logarithmic plot of the time dependence of the 
distance 1 — F from the target state |0) for the case where 
the system is initially in |00,0) and where k — 2r — 0.2 g 
and O = 0.03^. (b) Same as (a) but for a time dependent 
laser pulse with a Rabi frequency Q{t) as in Eq. (39) with 
Qo = 0.015^. 



the following that 



(l + 7c(0)t) 



2 ' 



(39) 



where 7c (0) denotes the cooling rate of the entangling 
scheme for the initial Rabi frequency Q{0) = CIq. Fig. 11 
confirms that choosing a time-dependent Rabi frequency 
Q indeed improves the speed of the entangling scheme 
without sacrificing much of its quality. 



IV. CONCLUSIONS 

In this paper, we propose an entangling scheme for two 
atoms trapped inside an optical cavity. Each atom should 
contain a A-like level configuration with the ground states 
|0) and |1) forming one qubit and an excited state |2) 
(cf. Fig. [T]). Three laser fields should be applied simul- 
taneously. Two of them continuously drive the 0-2 tran- 
sition which is in resonance with the cavity mode, while 
the third laser drives the 1-2 transition. Individual laser 
addressing of the atoms is not required. Most impor- 
tantly, the laser detunings should be chosen as proposed 
in Eq. (32) in Section HID As a result, the maximally 
entangled atomic ground state |+,0) with no photons 
in the cavity becomes the stationary state of the atom- 
cavity system. To complete the state preparation, the 
laser fields should be turned off after a certain transi- 
tion time. The presence of non-zero spontaneous decay 
channels, i.e. the leakage of photons through the cavity 
mirrors and direct spontaneous emission from the atoms, 
are essential for the scheme to work. 

The proposed state preparation scheme is a concrete 



realisation of a recent proposal to cool atoms into entan- 
gled sta te [23] . Choosing the laser detunings as proposed 
in Eq. (32) guarantees that only the target state |+,0) 
experiences off-resonant driving. All other states with 
no population in |2) and in the cavity mode interact res- 
onantly with one of the three applied laser fields. As 
in laser sideband cooling [29 , this makes it much more 
likely for the atom-cavity system to decay into the target 
state than being driven out of it. As a result, most of 
the population of the system accumulates in |+,0) with 
both atoms in a well-defined, maximally entangled state. 
Since the relevant detuning of this state is essentially 
given by the atom-cavity coupling constant ^, the anal- 
ogy to laser sideband cooling suggests that the scheme 
works best when all other system parameters are much 
smaller than g. 

Due to laser driving with three different laser fields, it 
is not possible to solve the time evolution of the proposed 
entangling scheme analytically. The reason is that there 
is no interaction picture in which the Hamiltonian of the 
system becomes time independent. To obtain at least 
approximate analytical solutions for the cooling rate 7c 
and the stationary state fidelity F (cf. Eq. ([37|)), Section 
|n| discusses a closely related state preparation scheme for 
a much simpler analytically tractable toy-model. A com- 
parison with this toy-model provides much insight into 
the state preparation via cooling as well as analytical re- 
sults. These are confirmed in Section IIIE| by extensive 
numerical solutions of the time evolution of the atom- 
cavity system in Fig. [l] Section IIIF| finally suggests a 
method to speed up the state preparation without sac- 
rificing its fidelity by using time-dependent laser fields 
with rapidly decreasing Rabi frequencies (cf. Eq. ([39|). 

Compared to other recent entangling schemes for 
atom-cavity systems [271 [2E| , the scheme proposed here 
predicts higher fidelities for the same experimental pa- 
rameters. As illustrated in Fig. 10 in Section [ill E[ it can 
achieve fidelities above 90% even when the single-atom 
cooperativity parameter C is as low as 20. With a co- 
operativity parameter C = 25 we can achieve a fidelity 
of 93%, while Ref. ^27j predicts fidelities above 92% only 
for C > 50. Ref. [27] uses a similar level scheme as our 
proposal but with the addition of a driven microwave 
transition between the triplet states. Ref. [28 requires 
the presence of a magnetic field gradient to produce the 
required level splittings to cool atoms into an entangled 
state. Compared to other quantum computing schemes 
using dissipation [7l-[9l |Tl] IT3HT7] , the state preparation 
scheme discussed here no longer relies on the detection 
of single photons or macroscopic fluorescence signals to 
herald the success of the state preparation. Its implemen- 
tation might already be in reach with current technology. 
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